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Abstract. In this paper, we use the approximation of shallow water waves 
(Margaritondo 2005 Eur. J. Phys. 26 401) to understand the behavior of a tsunami 
in a variable depth. We deduce the shallow water wave equation and the continuity 
equation that must be satisfied when a wave encounters a discontinuity in the sea 
depth. A short explanation about how the tsunami hit the west coast of India is given 
based on the refraction phenomenon. Our procedure also includes a simple numerical 
calculation suitable for undergraduate students in physics and engineering. 
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1. Introduction 

Tsunamis are water waves with long wavelengths that can be triggered by submarine 
earthquakes, landslides, volcanic eruption and large asteroid impacts. These non- 
dispersive waves can travel for thousands of kilometers from the disturbance area where 
they have been created with a minimum loss of energy. As any wave, tsunamis can be 
reflected, transmitted, refracted and diffracted. 

The physics of a tsunami can be very complex, especially if we consider its creation 
and behavior next to the beach, where it can break. However, since tsunamis are 
composed of waves with very large wavelengths, sometimes greater than 100 km, they 
can be considered as shallow waves, even in oceans with depths of a few kilometers. 

The shallow water approximation simplifies considerably the problem and still 
allows us to understand a lot of the physics of a tsunami. Using such approximation, 
Margaritondo pQ deduced the dispersion relation of tsunami waves extending a model 
developed by Behroozi and Podolefsky [2]. Since energy losses due to viscosity and 
friction at the bottom [3] can be neglected in the case of shallow waves, Margaritondo, 
considering energy conservation, explained the increase of the wave height when a 
tsunami approaches the coast, where the depth of the sea and the wave velocity are 
both reduced. 

In this paper we use one of the results of Ref. [I] in order to deduce the wave 
equation and include the variation of the seabed. Thus, we are able to explain the 
increase of the wave amplitude when passing from deeper to shallow water. Also, we 
discuss the refraction of tsunami waves. This phenomenon allowed the tsunami of 
December 24, 2004, created at the Bay of Bengal, to hit the west coast of India (a 
detailed description is given in 01)- These both inclusions - the seabed topography and 
the wave refraction - where pointed out by Chu |5. as necessary to understand some 
other phenomena observed in tsunamis. 

This paper is organized as follows. The wave equation and the water flux 
conservation are used in section 2 in order to explain how and how much a shallow 
wave increases when passing from a deeper to a shallow water. In section 3, we extend 
the results obtained in section 2 to study how a wave packet propagates in a water tank 
where the depth varies; in this section we use some numerical procedures that can be 
extended to the study of any wave propagating in a non-homogeneous medium. Also, 
the refraction of the 2004 tsunami in the south on India is discussed in section 3. The 
shallow wave and the continuity equations are deduced in appendix A. 

2. Reflection and transmission of waves in one dimension 

Consider a perturbation on the water surface in a rectangular tank with a constant 
depth. In the limit of large wavelengths and a small amplitude compared with the 
depth of the tank, the wave equation can be simplified to (see Appendix A) 




gh 




(1) 



dt 2 



dx 2 ' 
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where y(x, t) is the vertical displacement of the water surface at a time t, propagating 
in the x direction, g is the gravity acceleration and h is the water depth. 

Equation is the most common one- dimensional wave equation. It is a second- 
order linear partial differential equation and, since g and h are constants, any function 
y = f(x ± vt) is a solution (v = \fgh is the wave velocity). An interesting aspect of eq. 
(JTJ is that a propagating pulse does not present a dispersion due to the same velocity 
of all wavelengths and, thus, preserves its shape. Light in vacuum (and, in a very good 
approximation, in air) is non- dispersive. Also, sound waves in the air are nearly non- 
dispersive. (If dispersion was important in the propagation of the sound in the air, a 
sound would be heard different in different positions, i.e., music and conversation would 
be impossible) 

However, the velocity of shallow-water wave varies with the depth. Thus, shallow- 
water waves are dispersive in a non-uniform seadepth. 

In order to study the evolution of a tsunami in a rectangular box with variable 
depth, which will be detailed in the next section, we approximate the irregular depth by 
successive steps. So, in the next paragraphs we will explain the treatment used when a 
wave encounters a discontinuity. 

Every time the tsunami encounters a step, part is transmitted and part is reflected. 
Then, consider a wave with an amplitude given by y — cos(kx — out), where k and u are, 
respectively, the wave number and the frequency incoming in a region where the depth 
of the water, and also the wave velocity, have a discontinuity as represented in Fig. [T] 
On the left-side of the discontinuity the perturbation is given by 
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Figure 1. Evolution of a pulse propagating from a deep sea to a shallow sea (solid 
line). The dashed line is the wave velocity in units of 1/200 m/s. The upper frame (A) 
shows the pulse before the velocity discontinuity, the middle frame (B) at the velocity 
discontinuity and the lower frame (C) after the discontinuity. 



It 




2/i (x, t) = cos(kx — ujt) + Rcos(kx + ut + (fx), (2) 

where Rcos(kx + ut + <fi) corresponds to the reflected wave and <pi is a phase to be 
determined by the boundary conditions. On the right-side of the discontinuity the wave 
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amplitude is given by 

y 2 (x, t) = T cos(k'x — ut + (p 2 ), (3) 

corresponding to the transmitted wave part. The wave numbers for x < and x > 
are, respectively, 

* = £ (4, 

and 

k' = -„ (5) 
v' 

where v and v' are the velocities of the wavepacket at the left and right sides of the 
discontinuity. 

In order to determine R and T we must impose the boundary conditions at x — 0. 
For any instant, the wave should be continuous at x = 0: cosut + Rcos(ut + (pi) = 
Tcos(— ut + <f '2). The same should happen with the flux, f(x,t), given by (see eq. 

(EH) 

dz{x,t) 

f(x,t) = h — , (6) 

where z(x, £) is the horizontal displacement of a transversal section of water (see equation 
(jA.3j) for the relation between z and ?/). 

Imposing the boundary conditions yi(0,t) = 2/2(0, t) and /i(0, £) = /2(0, t) we can 
deduce smipi = sin^ 2 = 0. Then choosing tpx = <p 2 = we obtain 

fl = k '~ k = v ~ v ' , 7) 
k + k' v + v' 

and 

k + k' v + v' 

It is worthwhile to mention here that other choices of </?i and ^2 will change the 
signs of i? and T. However, in this case, it will also change the phases of the reflected 
and transmitted waves. Both modifications will compensate themselves and the shape 
of the wave will remain unchanged in relation of the choice (pi = f>2 = 0. 

Reflection and transmission are very important effects in wave propagation: every 
time a traveling wave (light, water waves, pulses in strings, etc.) encounters a 
discontinuity in the medium where it propagates, reflection and transmission occur. 
Since there is no energy losses, energy flux is conserved. The energy of a wave is 
proportional to the square of its amplitude [l] . Thus, the energy flux is proportional to 
the squared amplitude times the wave velocity. The energy flux of the incident wave at 
x = is given by 

(Pine = V (9) 

(the amplitude of the incident wave was chosen as 1). 

The reflected and transmitted energy flux are given by 

(f>refl = R 2 V (10) 
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and 



trans 




(11) 



respectively. 

Using eqs. (JJJ) and (JBJ) it is immediate to show that 

4*inc *firefl "I - *fitrans m 



(12) 



It is worthwhile to mention here that eqs. ( J 7 ]l . (|Sjl and (fT2*|) are classical textbook results 
on wave at interfaces. 

Fig. Q shows the evolution of a wavepacket in three distinct situations: before find 
a step, passing the step and after passing the step. On the left side of x = the depth is 
2000 m and on the right side 10 m. We can note a growth of the wavepacket amplitude 
after passing a step due to the velocity variation, in this case k'/k = 14, and, then 



Using energy conservation, Margaritondo deduced that the wave amplitude, when 
going from a sea depth hi to h 2 , increases by the factor [h\jh^) x ^ [1 . According to 
this result, the amplitude in our example should grow by a factor of 3.8. However, the 
growth observed was 1.9. The difference is due to the fact that when a wave packet 
encounters a discrete step of velocity, part of the energy is reflected (the reader can 
verify that eq. ([12)) is satisfied). As will be shown in the next section, when the sea 
depth varies smoothly, the reflected wave can be neglected, and our results becomes 
equal to Margaritondo's result. 

3. Waves in a variable depth 

In order to study the evolution of a tsunami when it propagates in a rectangular box 
where the depth varies, we initially made a wavepacket propagating into a crescent x 
direction. The variable depth was approximated by a succession of steps taken as narrow 
as we wish. 

The evolution of the wave packet was calculated as follows: 

• At time t, the wave packet amplitude y(x,t) was divided in n small discrete 
transversal sections of length Ax (in our case Ax = 50000 m). 

• Every small part of the wave packet y(x,t) was investigated: 

(i) If in a time interval At it stays in the same velocity step, then the wave 
packet at t + At was simply increased by y(x, t) at the position x + v(x)At. 

(ii) If in a time interval At, we choose At as 30 s, it encounters a velocity step, 
part is reflected and part is transmitted. The reflected and transmitted parts were 
calculated from eqs. (jZj) and (JBJ and added to the wave packet at t + At propagating 
to the left or to the right, respectively. The step width and the time interval At 
were chosen such that never the reflected or transmitted parts encounter a second 
step. 



T 



1.9. 
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Figure 2. The upper frame (A) shows an initial wave packet traveling to the right. 
The frames (B) and (C) show, respectively, the tsunami in an intermediary position 
and near to the coast. The tsunami velocity is given by the dashed line. Note that 
the wavepacket extension diminishes in the same proportion as the velocity. (Only the 
progressive part of the wavepacket is shown) . 



Fig. El shows three positions of the right-propagating wavepacket (the left- 
propagating was omitted). The initial wavepacket, fig. 2A, has its center at a depth 
of about 3930 m (v ~ 200 m/s), then the center of the wavepacket goes to a position 
where the depth is about 165 m (v ~ 40 m/s), fig. 2C. The growth of the amplitude 
is about 1.7. The difference between our result and the one expected by Ref. pQ 
((3903/165) 025 = 2.2) is due to the fact that we approximate the continuous depth 
variation by discretes steps and, as a consequence, the left-propagating wavepacket was 
not negligible. 

In the last paragraphs of this section we insert a short discussion of the refraction 
phenomenon. 
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Figure 3. Change of orientation of the wave crests (solid lines) in a variable sea depth. 
The dashed lines are the trajectories of the front waves. The sea depth varies with the 
q coordinate. 

Consider, for instance, a water wave propagating in a medium where the sea depth 
varies with q, as shown in fig. |3] (for instance q can be the distance from the coast). 
The wave crest at qi has a velocity v\ and at q2 velocity It is a matter of geometry 
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to show that a wavefront will change orientation as making a curve with radius 



dq 

For instance, consider what happened in south of India. The seadepth varies from 
about 2000 m at q\ =500 km far from the coast to about 100 m near the coast, q<i ~ 0. 
Thus, \dv/dq\ w 2.2 x 10~ 4 s -1 . As a consequence, the radius formed by a wave crest 
varies from about 640 km far from the coast to about 140 km near the coast. This is 
about what we can see from fig. 4. In fig. 01 a tsunami wavefront propagation, part in 
deep water and part in shallow water near the coast, will refract and, in consequence, 
change orientation. Fig. 0] shows a refraction map for the December 26, 2005 tsunami. 
The dashed curves are the frontwaves at 100, 150, 200 and 300 minutes after the 
earthquake [I]. 




Figure 4. The dashed curves are the wavepackets corresponding, approximately, to 
the tsunami of 26 December 2004. The star shows the earthquake epicenter. 



4. Discussion 

The tsunami of 26 December 2004, obviously, did not propagate in a rectangular 
box, but, approximately, in a circular tank. As the tsunami propagates, its extension 
increases and, in consequence, its amplitude diminishes. However, when it approaches 
shallow waters, near to the coast, its amplitude grew up again as shown by the simplified 
model developed in this paper. 

The developed model depends on two approximations: the wave amplitude is small 
and the length of the wavepacket is large when compared with the sea depth. Since the 
first approximation is not valid near the beach, we stopped the evolution when the front 
part of the tsunami wavepacket attained a depth of about 50 m, shown in Fig. |2P- 

In summary, with the model developed in this paper we showed how to use the 
approximation of shallow water waves to a variable depth. This simple model allows 
us to understand what occurs when a tsunami goes from deeper to shallow waters: the 
velocity of the rear part of the wavepacket is larger than the velocity of its front part, 
causing water to pile up. Also, refraction effects, that are not present in a sea of constant 
depth, can be observed near to the coast. 



Understanding the tsunami with a simple model 



8 



MTY thanks the Brazilian agency FAPESP (Fundagao de Amparo a Pesquisa do 
Estado de Sao Paulo) for financial support. 

Appendix A. Deduction of the wave equation for waves with wavelengths 
much greater than the water depth 

We will deduce the wave equation in a simplified situation. We will make the following 
approximations (all of them can be applied to tsunamis located far from the beach): the 
part of the restoration force that depends on the surface tension can be neglected in the 
case of waves with large wavelengths; the wavelength or the extension of the wavepacket 
will be considered much longer than the depth of the water (in the case of tsunamis the 
wavelengths and the ocean depth can have, approximately, hundreds of km and a few 
km, respectively); the wave amplitude will be considered much smaller than the ocean 
depth. Another simplification is the tank where the wave propagates: we will consider 
a wave propagating in a rectangular box with vertical walls and constant depth. In this 
approximation of shallow water waves, all the droplets in the same transversal portion 
have the same oscillatory horizontal motion along the x direction. Finally, friction at 
the bottom will be neglected 

Fig. IA1I illustrates the situation considered. The wave direction of propagation is 
x; h Q is the unperturbed height of the water; L is the box width. 
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Figure Al. Rectangular box with vertical walls and constant depth where the wave 
propagates. 
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Figure A2. Side-view of Fig. IA1I showing a perturbation in the water. 



Fig. IA2I shows the same box of Fig. IA1I in a side-view showing a perturbation in 
the water. A lamellar slice with width Ax in x and a height h at a time t, will have a 
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height h + y and a width Ax + Az when it occupies the position x + z at an instant 
t + At. z = z(x, t) is the horizontal displacement - along the x direction - of a vertical 
lamellar slice with an equilibrium position at x. When the wave propagates, this part of 
the water oscillates to left and right. Equaling the volume of water in Ax and Ax + Az, 
we have: 

LhAx = L(h + y){Ax + Az) 

= L(hAx + hAz + yAx + yAz). (A.l) 

If we consider y « h and Az « Ax, then eq. (jA.l|) becomes 

hAz + yAx = 0, (A.2) 

or 

„=-k£. (A.3) 

This last equation is the mass conservation equation of the fluid and relates the vertical 
displacement of the water surface, y, with the horizontal displacement of a vertical slice 
of water. 

To apply the second Newton Law to a small portion of the lamellar slice, Ah, of 
water (see Fig. IA2J) . we should calculate the total force, AF, acting on it. This force 
depends on the pressure difference between the opposite sides of the slice: 

dP 

AF = AhL (P(x) - Fix + Ax)) ~ -AhL—Ax. (A.4) 

ox 

Then, F = ma leads to 

dP d 2 z 
- AhL—Ax = pLAhAx—, (A. 5) 

ox ot 2 

where m is the mass of the water slice, a is the acceleration of the transversal section 

of water given by the second partial derivative of z with respect to t, and p is the water 

density. 

Since 

dP dy 

dx ^ dx : ^ ^ 

where g is the gravity acceleration, we can write eq. (|A.5J) as 

d 2 z dy 

W = ~ g Tx (A ' 7) 

Derivate both sides of eq. (jA.3|) with respect to x we have 



dx ^ dx 2 ^ ^ 

Finally, using eq. (jA.8|) in eq. (jA.7|) we obtain the wave equation [HJ |E] 

d 2 z , d 2 z , . s 

gh— (A.9) 



dt 2 dx 
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Using eq. (jA.3|) we can show that the vertical displacement of the water surface, y, 
obeys an equivalent wave equation: 

<-°> 

The solutions of eqs. (|A.9J) and (|A.10|) are any function of x — v t or x + vt, where 
the wave velocity is given by 

v = \fgh. (A.ll) 

Eq. (jA.llJl is a particular case of the general expression for the dispersion relation 
for waves in water surface 



/ qX 27rcr\ , 2nh . . . 

\{h + ^) Unh — (A12) 

where A is the wavelength, p the water density and a is the water surface tension. In 
the case of long wavelength and neglecting surface tension, eq. (jA.12j) reduces to eq. 

4HU). 

Eq. (jA.ll|) gives two useful conclusions for waves presenting the same characteristics 
of a tsunami: the wave velocity does not depend on the wavelength, and the wavepacket 
does not disperse when propagating in a region of constant depth. 

Since z = z(x, t) is the horizontal displacement of a lamellar slice of water, then 
the water flux, /, is given by 

/ = Lh^. (A.13) 



References 

[1] Margaritondo G 2005 Eur. J. Phys. 26 401 

[2] Bchroozi F and Podolefsky N 2001 Eur. J. Phys. 23 225 

[3] Crawford F S, 1987 Am. J. Phys. 55 171; Waves - Berkeley physics course volume 3 (McGraw-Hill 
book company, 1968). 

[4] Lay T et al., "The great Sumatra-Andaman earthquake of 26 December 2004", Science 308, 1127- 

1133 (2005). 
[5] Chu A K H, Eur. J. Phys. 26 L19 

[6] Alonso M and Finn E J, Physics (Addison Wesley, 1992). 



